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1 Abstract 

The effect of the anisotropy in the exchange interaction and film thickness on the hysteresis behavior 
of the anisotropic Heisenberg thin film has been investigated with effective field formulation in a 
two spin cluster using the decoupling approximation. The behaviors of the hysteresis loop area, 
coercive field and remanent magnetization with the film thickness and anisotropy in the exchange 
interaction have been obtained. 
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2 Introduction 

The magnetic properties of thin films are drastically different from the bulk counterparts. In the 
presence of the free surfaces the translational symmetry of the system broken, due to the surface 
atoms are embedded in an environment of lower symmetry than that of the inner atoms mm- 
The most important result of this fact that, the surface of the thin film could show the ordered 
phase even if the bulk is in the disordered phase. Depending on the value of the surface exchange 
interaction, this interesting behavior may occur and this fact has been observed experimentally [31 
HE]- Another interesting magnetic behavior of thin films is dependence of the critical temperature 
and saturation magnetization value with film thickness. It was shown that the critical temperature 
and the average magnetic moment per atom increases with the increasing thickness of the film HE]- 

One class of the films which exhibit a strong uniaxial anisotropy [8] can be modeled by Ising 
model, which is the highly anisotropic limit of the Heisenberg model. This type of model on thin 
film geometry has been widely studied in literature by means of several theoretical methods. Widely 
used methods for these systems are Monte Carlo (MC) simulations [9], mean field approximation 
(MFA) |T[)] and effective field theory (EFT) ITT] . In order to mimic the surface effects on the 
system, different exchange interactions have been defined and the model has been solved within 
the framework of EFT [12] . Also more realistic effects such as amorphisation of the surface due to 
environment have been handled within the same framework m- There are also higher spin Ising 
thin films, e.g. spin-1 Ising thin films have been studied [141 115] . 

On the other hand, wide variety of materials do not have strong uniaxial anisotropy. In this 
case, the model has to include the spin-spin interactions as Heisenberg model. The most simple 
system with surfaces is semi infinite model, which has only one surface, in contrast to the thin film. 
Heisenberg model on a semi infinite geometry has been solved using a wide variety of techniques 
such as Green function method [16] , renormalization group technique [17], MFA [18] , EFT [T9j [20] 
[21] . high temperature series expansion [221123] , Besides, critical and thermodynamic properties of 
the bilayer [24] [25] and multilayer [26] systems have been investigated within the cluster variational 
method in the pair approximation. Heisenberg model in a thin film geometry has been solved in 
a limited case. For instance, Green function method [27] 2s, [26] . renormalization group technique 
[30 j, EFT [HU [32] and MC [33l [34] , are among them. 

Hysteresis is a common behavior of the most of the physical systems. It originates from the delay 
of the response of the system to the driving cyclic force, and shows itself as a history dependent 
response. Magnetic hysteresis is one of the most important and interesting features of the magnetic 
materials. Hysteresis loop area (HLA), coercive field (CF) and remanent magnetization (RM) are 
the parameters that give important clues about the shape of the hysteresis loops. Also, these 
hysteresis properties are very important in technological applications such as manufacturing of 
magnetic recording media. The RM is defined as residual magnetization which is the magnetization 
left behind in the system after an external magnetic field is removed. CF is defined as the value 
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of the external magnetic field needed to reverse the sign of the magnetization. On the other 
hand, HLA is simply the area of the hysteresis loop, which corresponds to energy loss due to the 
hysteresis. In our recent work we obtain to the hysteresis behaviors of the anisotropic Heisenberg 
model on the bulk system within the EFT formulation [35]. 

The aim of this work is the determine the effect of the anisotropy in the exchange interaction 
and number of layers of the thin film, on the hysteresis behavior of the Heisenberg thin film. The 
phase diagrams and the magnetization behaviors of this model has already been obtained [551 . 

For this aim, the paper is organized as follows: In Sec. [3] we briefly present the model and 
formulation. The results and discussions are presented in Sec. [2 and finally Sec. [5] contains our 
conclusions. 


3 Model and Formulation 
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Figure 1: Schematic representation of the thin film. 

The schematic representation of the thin film can be seen in Fig. [l] System is infinitely long in 
x and y directions, while finite in z direction. Thin film can be treated as a layered structure which 
consist of interacting L parallel layers. Each layer (in xy plane) is defined as a regular lattice with 
coordination number z. When we choose z = 4, this means that each layer is a square lattice and 
each nearest neighbor layer have interaction. The Hamiltonian of the thin film is given by 

n = - E [ A d + s\s)) + st*3] - (!) 

<i,j> i 

where s£,sf and sf denote the Pauli spin operators at a site i. Jij is the exchange interaction 
and Aij stands for the anisotropy in the exchange interactions between the nearest neighbor spins 
located at sites i and j. Longitudinal magnetic field at a site i is denoted by Hi in Eq. (JTJ) . The 
first sum is carried over the nearest neighbors of the lattice, while the second one is over all the 
lattice sites. The exchange interaction and anisotropy in exchange interaction between the spins 
on the sites i and j take the values according to the positions of the nearest neighbor spins. The 
two surfaces of the film have the intralayer interactions (Ji, Ai). The interlayer coupling between 
the surface and its adjacent layer (i.e. layers 1, 2 and L — 1,L) is denoted by (J 2 , A 2 ). For the rest 
of the layers, the interlayer and the intralayer couplings are assumed as (J 2 , A 2 ). 

In order to get the hysteresis behavior of the anisotropic Heisenberg thin film, we use the EFT, 
which can provide results that are superior to those obtained within the MFA. This fact comes from 
the consideration of self spin correlations, which are omitted in the MFA. In an EFT approximation, 
one constructs a finite cluster. The anisotropic character of the spin-spin interactions cannot be 
handled in one spin cluster, thus we use two spin cluster approximation here. EFT for the two 
spin cluster namely EFT-2 formulation m , was first proposed in Ref. [38] for Ising systems. In 
EFT-2 approximation, we choose two spins (namely si and s 2 ) in each layer and treat interactions 
exactly in this two spin cluster. In order to avoid some mathematical difficulties, we replace the 
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perimeter spins of the two spin cluster by Ising spins (axial approximation) |39j . After all, by 
using the differential operator technique and decoupling approximation (DA) [5D], we can get an 
expression for the magnetization per spin, i.e. m = (sf + s l))- I n the thin film geometry, the 
number of L different representative magnetizations for the system (by following the procedure 
given in Ref. E21) can be given as, 


m l = (©l,l® 2 > 2 )/l ( x ^y^Hl,H 2 ) |x=0,3/=0 

m k = ^0 2 ,fe—i02,fe©2,fc+i^ fi (x,y,H 1 ,H 2 ) \ x =o, y =Oik = 2,3,..., L - 1 (2) 

m L = (02,l-iB? iL ) fi (x,y,H 1 ,H 2 ) \ x =o, y =o- 

Here m*, (i = 1, 2,..., L) denotes the magnetization of the i th layer. The operators in Eq. © are 
defined via 

— [Afca, “b TTllBkx\ [Afcy + Pk^lBky] (3) 

where 

B km = sinh (J|V m ), k = 1,2; m = x, y. ^ ' 

The functions in Eq. © are given by 


fn (x,y,Hx,H 2 ) 


x + y + Hi + H 2 

X^ n) cosh 


sinh 

+ exp (-2 PJ%) cosh (/3F 0 (n) ) 


where 


X, 


(n) 


Y r 


( n ) 


(x + y + H 1 + H 2 ) 

(2A„ J^) 2 + (x — y + Hi — H 2 ) 2 


(5) 

( 6 ) 


with the values n— 1,2. In Eq. ©> we set /3 = 1 /(fc^T) where kg is Boltzmann constant and T 
is the temperature. 

Magnetization expressions given in closed form in Eq. © can be constructed via acting differ¬ 
ential operators on related functions. 

Differential operator technique depends on writing Eq. © in a polynomial form in m via the 
exponential differential operator. The effect of the exponential differential operator to an arbitrary 
function F{x) is given by 

exp (aV)F (x) = F (x + a) (7) 


with any constant a. 

By using Eq. © in Eq. © and converting the hypertrigonometric functions to exponentials, 
after using the Binomial expansion, Eq. © can be written in the form 
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where the coefficients are defined by 
Ki {p, q, r, s ) 


/l3-p d-9/l-r- .l-s, 
^T -1 „ -Alo 
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Mcc ^lp '^ 1 2rc "^2 p 


K 2 (p, q, r, s, t, v) = \ r J y s J ^ 2 * 

2 p+r+t u q+s+v 


a: 


5 —(p+r+t) A-iq+s+v) , 


^1/ 


B 2x B%y S+v f 2 (x, y, Hi,H 2 ) |x=o,»=o. 


( 8 ) 


(9) 


3 





These coefficients can be calculated from definitions given in Eq. © by using Eqs. © and ©• 
Eq. © is a system of coupled non linear equation system, and it can be solved via usual 
Newton-Raphson iteration [41]. The solution of the system is the longitudinal magnetizations of 
each layer ( rrii,i = 1,2 ,,L). The total longitudinal magnetization (m) can be calculated via 



i—1 


( 10 ) 


We can determine the second order critical point by linearizing the equation system given in 
Eq. ©. Since all the longitudinal magnetizations are close to zero in the vicinity of the second 
order critical point, the solution of the 
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linear equation system will give the second order critical point. 


4 Results and Discussion 


The parameters (Jf, Ai) give the exchange interaction of the nearest neighbor spins at the surfaces, 
while (J|, A 2 ) gives the remaining exchange interactions between the nearest neighbor spins of the 
thin film. 

Let us select the unit of energy as J (J > 0) and scale the exchange interaction with J as, 


Ti = 


J[_ 

J ’ 


* = 1,2 


( 12 ) 


This defined dimensionless parameters in Eq. oa are all positive or zero, since the system has 
only ferromagnetic interactions. Let the parameter q to control the ratio of the anisotropy in the 
exchange interactions between the surfaces and other layers, 


Q = 


Ai 

A2 


(13) 


4.1 Isotropic Model 

In this case, 

Ai = A 2 = 1.0. (14) 

In order to see the effect of the ratio of r\/r 2 on the hysteresis behavior of the film, let us choose 
r 2 = 1.0. The phase diagram of this system has already been obtained [3B]. The phase diagram 
of the thin film that have isotropical exchange interaction shows well known interesting behavior. 
The phase diagrams for different film thickness ( L ) intersect at a special point which can be 
denoted by (r^ksT*/J) in the (ri,fcsT c /J) plane. This fact comes from the surface of the 
system. Magnetically disordered surface can coexist with a magnetically ordered bulk phase for 
the values of ?q that provide rq < r\ while for the values rq > , surface can reach the magnetically 

ordered phase before the bulk. The special point coordinate of the thin film has been obtained as 
(7-*, ksT*/J) = (1.345,4.891) within the EFT-2 formulation [36], which is nothing but the values 
that makes the critical temperature of the film independent of the thickness of the film L. This 
critical temperature value is just the critical temperature of the corresponding bulk system (the 
system with simple cubic lattice) within the same model m- 

In order to see the effect of the film thickness L on the hysteresis behavior, we choose four 
different values of the (7q, fc^T/J), while r 2 = 1.0 as shown in Fig. [2] These behaviors can be seen 
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Figure 2: The selected values of the temperatures and rq for the inspection of the hysteresis 
behaviors in the isotropic model, shown in the phase diagrams of the thin film in isotropic case. 


in Fig. [3] for film thicknesses L = 3,4,6,10. Typical behavior for the hysteresis loops with rising 
thickness for n < r* is as expected and can be seen in Figs. [3] (a) and (b). When the film gets 
thicker, ferromagnetic hysteresis loops appear, since rising thickness induce a ferromagnetic phase, 
as seen in Fig. [2] In a ferromagnetic region, rising thickness give rise to wider hysteresis loops (e.g. 
compare curves labelled by C and D in Fig. [3] (a)). The behavior of the RM, CF and HLA can 
be seen more clearly in Fig. [4] for this region. The parameter values in Fig. [4] are rq = 0.2 < r*, 
rq = 1.0 and /c^T/J = 3.5. For this set of values, only the film that has thickness L = 3 is in the 
paramagnetic phase, as seen in Fig. [2] Rising thickness, first rises the RM, CF and HLA smoothly. 
Then, there is no significant change of these values, when the film gets thicker. 

At a special point (r^ksT* / J) = (1.345,4.891), rising film thickness does not change the 
hysteresis loops. At this value, all hysteresis loops are the same, as seen in Fig. [3] (c). Lastly, 
when the parameter values selected as (rq, ksT/J) = (1.8, 5.3), rising film thickness effects on the 
hysteresis loops reversely in comparison with the parameters that lies on the left of the special 
point. As seen in Fig. [3] (d), hysteresis loops get narrower when the film thickness rises. 

All these behaviors can be explained by the special behavior of the thin films namely ordinary 
and extraordinary phase transition behaviors. For the values that rq < bulk part of the film 
is dominant for the magnetic properties of the film, where ordinary phase transition takes place 
and thicker films have higher critical temperatures. On the other hand for the values that provide 
rq > r\ 1 extraordinary transition takes place. In this case surface is dominant, and thicker films 
have lower critical temperatures and the surface of the film can be ordered before the bulk when 
the temperature is lowered. Also at any temperature which is below the critical temperature and 
not so close to the zero, the value of the surface magnetization is reduced in comparison with the 
bulk for the rq < r^, and vice versa. This facts are valid for the isotropic Heisenberg model [36] 
as well as the Ising model [12] and this outcome was also proven by more sophisticated techniques 
such as series expansion methods and Monte Carlo simulations [42] . We can conclude for this 
section that, in the ordinary region (i.e. the parameter values that the system can display ordinary 
transition, rq < r^) rising film thickness enhances the CF, RM and HLA, while the reverse is true 
for the extraordinary region (i.e. rq > r^). 

4.2 Anisotropic Model 

In this case the anisotropies in the exchange interactions are related to each other as given in Eq. 
da and A 2 = 1.0, i.e. the interior of the film consists of completely isotropic exchange interactions 
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Figure 3: Hysteresis loops for the isotroic Heisenberg thin film for the r 2 = 1.0 and thicknesses L = 
3,4,6,10, for selected values of the (n, ksT/ J) = (0.5,4.2), (1.2,4.6), (1.345,4.891) and (1.8, 5.3). 


while the surfaces of the film have the exchange interactions of the Ising type anisotropic {q = 0), 
completely isotropic ( q = 1) or XXZ type anisotropic (0 < q < 1), according to the value of q. 

Hysteresis loops of the anisotropic thin film can be seen in Fig. [5] While Fig. 0(a) corresponds 
to the Ising type surface and Fig. 0(d) isotropic Heisenberg type surface, Figs. 0(b) and (c) are 
related to the XXZ type anisotropic surface with the values of q = 0.3 and q = 0.7, respectively. 
The temperature is chosen as fcgT/ J = 4.3 in Fig. 0 

First, all figures exhibit the same behavior with rising film thickness, that is, rising film thickness 
enlarges the hysteresis loops. This regular behavior is valid for the cases with different types 
of surfaces (i.e. surfaces that have different anisotropy in the exchange interaction). The only 
difference is, the loop labeled by A in Fig. 0 (d) which corresponds to the paramagnetic phase for 
the film thickness L = 3, while all other loops are related to the ferromagnetic phase. 

On the other hand, when the anisotropy in the exchange interaction of the surfaces of the film 
is lowered, then the hysteresis loops become narrower (e.g. see the loops labeled by D in Figs. 0 
(a), (b), (c) and (d), which are the loops for the film thickness L = 6 for different anisotropy in the 
exchange interaction of the surfaces of the film). In other words, rising film thickness and rising 
anisotropy in the exchange interaction of the surfaces (i.e. decreasing q), affect the hysteresis loops 
in the same way, i.e. enlarging the loops. 

The shrinking behavior of the hysteresis loops with decreasing film thickness shows itself in the 
behavior of the HLA which can be seen in Fig. 0 (c), which is the variation of the HLA with q for 
selected values of film thickness L = 3,4,6,10 and temperature fc^T/J = 3.3. 

As shown in Fig. 0 (c), the regular decreasing behavior of the HLA with rising q turns a 
behavior of staying almost constant. The difference between the HLA of the films that have 
different thickness, is while the HLA of the film that have thickness L = 3 getting zero after a 
specific value of q (see curve labeled by A in Fig. 0(c)), the other curves cannot have the value of 
zero. 

The difference between the HLA values of the film that have different thickness is related to 
the CF as well as the RM. This fact can be seen in Figs. 0 (a) and (b). In Fig. 0 (a) we depict 
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Figure 4: Variation of the (a) RM, (b) CF and (c) HLA with the number of layers (film thickness) 
for selected values of ksT/J = 3.5, T 2 = 1.0 and r 2 = 0.2. 


the variation of the RM with q for selected values of L = 3,4, 6,10 and temperature ksT/ J = 3.3, 
while the same can be seen in Fig. [G] (b) for the CF. It can be seen from these Figs, that the 
thinner film can have paramagnetic phase when the anisotropy in the exchange interaction of the 
surface changes. 

5 Conclusion 

In this work, the effect of the anisotropy in the exchange interaction and the thickness of the film 
on the hysteresis behavior of thin films has been investigated. As a formulation, the differential 
operator technique and DA within the EFT-2 formulation has been used. 

The isotropic model on the thin film displays interesting behavior namely, for the values that 
provide r\ < rj when the thickness of the film rises, hysteresis loops get wider. This relation gets 
reversed for the values n > r*. At a special point (i.e. n = r\ and ksT/J = 4.891) the hysteresis 
loops are independent of the film thickness. 

On the other hand, in the anisotropic case the dependence of the hysteresis loops on the film 
thickness is more simple than the isotropic case. When the films get thicker, hysteresis loops get 
wider. This inspection has been made for several different cases by taking the interior of the film as 
isotropically interacting Heisenberg spins. The surface of the film that have spins have interaction 
of the Ising type, anistoropic Heisenberg type (by means of XXZ model) and isotropic Heisenberg 
type. One conclusion was, when the surfaces of the film have isotropically interacting Heisenberg 
spins, hysteresis loops are narrower than the other type of surfaces, when the other parameters 
(such as the temperature) are kept fixed. This is because of that, in the former case it is easier to 
follow the magnetic field than the latter case due to the isotropical spin-spin interaction. 

After q > 1.0 hysteresis loops of the thinner film (L = 3) disappear due to the transition from 
the ferromagnetic phase to the paramagnetic one, while the thicker films can stay in an ordered 
phase. 

We hope that the results obtained in this work may be beneficial form both theoretical and 
experimental point of view. 
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Figure 5: Hysteresis loops for the anisotroic Heisenberg thin film for the r 2 = 1.0, fc^T/J = 4.3 
and thicknesses L = 3,4, 6,10, for selected values of the (a) q = 0.0, (b) q = 0.3, (a) q = 0.7, and 
(a) q = 1.0. 
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